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Abstract 

Form factors of rare B i^Q (1430)£'''£~ decay are calculated within three-point 
QCD sum rules, with K^{1^30) being the p-wave scalar meson. The branching ra- 
tios are estimated when only short, as well as short and long distance effects, are 
taken into account. It is obtained that the B /Tq (1430)^^^^ {i = e,/i) decay is 
measurable at LHC. Measurement of these branching ratios for the semileptonic rare 
B — > 7^0(1430)^+^^ can give valuable information about the nature of scalar meson 
i^^(1430). 
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1 Introduction 



The flavor changing neutral current (FCNC) processes induced by & — > s{d) transitions 
provide the most sensitive and stringiest test for the standard model (SM) at one loop 
level, since they are forbidden in SM at tree level [1,2]. For this reason these decays are 
very sensitive to the physics beyond the SM via the influence of new particles in the loops. 

Despite the branching ratios of FCNC decays are small in the SM, quite intriguing 
results are obtained in ongoing experiments. The inclusive B Xsi~^i~ decay is observed 
in BaBaR [3] and Belle collaborations. These collaborations also measured exclusive modes 
B — > K£'^£~ [4-6] and B K*i~^£'' [7]. The experimental results on these decays are in a 
good agreement with theoretical estimations [8-10]. 

There is another class of rare decays induced by 6 ^ s transition, such as B ^ 
Kq2{^4:30)£^£^ in which B meson decays into p-wave scalar meson. The decays B — > 
K*{U30)£+£- and B K*{U30)£n- are studied in [11] and [12], respectively The main 
problem in these studies is the calculation of the transition form factors. Transition form 
factors of these decays in the framework of light front quark model [12] are estimated in 
[13] and [14], respectively. 

In the present work we calculate the transition form factors of B ^ Kq{1A30)£^£^ 
decay in 3-point QCD sum rules method [15] (for a review on the QCD sum rules method, 
sec [16]), where i^^^(1430) is the p-wave scalar meson. The main reason for studying 
B Kq{1A30)£~^£~ decay form factors in the framework of 3-point QCD sum rules method 
is that part of the form factors, namely, those appearing in calculation of the matrix clement 
of the axial vector current between initial and final meson states, can easily be obtained 
from the results for the Bg — > -Ds(2317) transition form factors [21], which are calculated 
within the same framework, with the help of appropriate replacements (see section 2). 
The quark structure of the scalar mesons have not been unambiguously determined yet 
and is still under discussion. In order to establish the inner structure of scalar mesons 
much more effort from theoretical and experimental sides are needed. This method is 
employed to calculate the form factors of heavy meson semileptonic decays, D — > K^eV^ 
[17], B D{D*)£Pe [18], D K{K*)ePe [19], D nePe [20], B K{K*)£+£- (first 
reference in [8]), etc. 

The work is organized as follows: In section 2, we calculate the relevant hadronic form 
factors of the B — > Kq{14:30)£~^£~ decay with the help of 3-point QCD sum rules. Section 
3 is devoted to the numerical analysis and discussion of the considered decay and our 
conclusions. 



2 Form factors of the B K^{im)i+i- decay 

The exclusive B Kq{1A30)£^£^ decay is described at quark level by 6 ^ s£^£^ transition. 
The effective Hamiltonian responsible for the b — > s£^£^ transition in the standard model 
is 



GFoVt,V* 

^eff - 
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TT 
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where C7, Cg'^^ and Cio are the Wilson coefficients. It is well known that the Wilson 
coefficient Cg^-^ has a perturbative part as well as a resonant part which comes from the 
long-distance effects due to the conversion of the real cc into lepton pair (see, for example 
the second reference in [8]) which can be written as 

The explicit expressions of C7, Cg^^ and Ciq can be found in [1] . Cg'^^ is usually parametrized 
by using Breit-Wigner ansatz, 

rres _ 37r (0) ^ ^ T{V, i+i-)mv^ 
^g ^ ~2 « 2 2 '■ ' 

« v,=V(il>-V(6.) "^v^i - ^ ^myTy, 

where a is the fine structure constant and C^^^ — 0.362. 

The phenomenological factors aSj for the B — > K{K*)i'^£~ decay can be determined 
from the condition that they should reproduce correct branching ratio relation 

B{B J/il^K{K*) K{K*)tr) = B{B J/ijK{K*))B{J/ij tr) , 

where the right-hand side is determined from experiments. Using the experimental values 
of the branching ratios for the B — > ViK{K*) and Vi — > £^£^ decays, for the lowest two 
J/ip and ip' resonances, the factor ae takes the values: aei = 2.7, ae2 = 3.51 (for K meson), 
and asi = 1.65, ae2 = 2.36 (for K* meson). The values of aSj used for higher resonances are 
usually the average of the values obtained for the J ftp and ip' resonances. But, unfortunately 
the mode B J/-?/'_ft'Q(1430) for the B Xq(1430)£+£" decay has not been seen yet and 
therefore the right-hand side of the above-mentioned equation for the B K^{lA30)i'^i~ 
decay is unknown. For this reason we are not able to determine the factors ae^ for the 
B K*{U30)i+i- decay. 

As we have already mentioned above, the values of the factors ae, are of order ae ~ 1 
or ae ~ 2, for the B — ^ Ki~^i~ and B K*£~'~i~ decays. For this reason, even though 
qualitatively, in order to determine the branching ratio for the B Kq{1A30)£^£^ decay 
with the inclusion of long distance effects, the values of ae^ we use are borrowed from 
B — > K£'^£~ and B — > K*£'^£~ decays, namely, we choose ae = 1 and ae = 2 and performed 
numerical calculations with these values. 

The amplitude for the B Kq{1430)£^£~ decay can be obtained after evaluating the 
matrix elements of the quark operators in Eq. (1) between initial B and final Kq{1A30) 
meson states. It follows from Eq. (1) that in order to obtain the amplitude for the B 
X^(1430)^+^- decay the matrix elements {K^ |s7^(l - 75)! S) and {K^ \siai,^q^{l + 75)] B) 
are needed. These matrix elements are parametrized in terms of the form factors as follows: 

(i^*(1430)(pO |s7m75&| B{p)) = U{q')r, + f.{q^)q, , (2) 
(X*(1430)(p') \sia^,q''^,h\ B{p)) = ^ [V,q' - (m| - m')q^] , (3) 
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where Vn 



p')^ and Qfj, = {p ~ p')n- For calculation of these form factors in the 
framework of the QCD sum rules method, we consider the following three-point correlators 



TJ-o(,^)j^(0)J^ 











(4) 
(5) 



where J^o = ds, — hi^^d and J^^ — sia^yh are the interpolating currents of the 
scalar Kq{1A2>Q), and B mesons and weak flavor changing quark currents, respectively. The 
expressions of the form factors f+{q^) and f-{q^) can be obtained from the results for the 
Bg Dsg(2317) transition form factors given in [21], with the help of the replacements 
Jd^^ fx*, ^DsQ tTT-K*, ^ mg, and — > — 0. For completeness we also present 
the sum rules for the form factors f+{q^) and f-{q^) below. 

For this reason wc proceed by describing in full detail the calculation of Ft by considering 
the correlator given in Eq. (5). This correlator can be decomposed into a set of the following 
independent Lorentz structures: 



(6) 



where 11 and n'-"^ are functions of p^, p'"^ and g^, and a|]J,^ are the other tensors which can 
be built using the vectors p and p' and the metric tensor g^y. 

In further analysis we will be interested only on the invariant amplitude 11^, for whose 
amplitude we write the dispersion relation 



IlT{.p\p'\q 



{2n) 



]^ poo pc 
n)^ J mi Jm 



ds' 



{s — p'^){s' — p'^) 



+ subtraction terms , 
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where s', Q"^) is the spectral density, Q"^ = —q^. According to the main idea of the 
QCD sum rules method, left hand side of Eq. (7) must be calculated at large Euclidean 
momenta p"^ and p'"^ with the help of operator product expansion (OPE). The right hand 
side of Eq. (7) is determined by saturating it with the lowest mesonic states. Applying 
then the Borcl transformation on the variables p^ and p'"^ suppresses the higher resonance 
contributions and higher power corrections, as a result of which we get the sum rules for 
the corresponding quantities. 

Let us start by calculating the hadronic part of the correlator (5). Saturating this 
correlator with B and K^{1A:2>Q) mesons and selecting the structure {p^p',, — PuP'^), we get 



Mq') 



rriB + rriK* rub {m\;» - p''^){m\ - p^) 



In deriving Eq. (8) we use 







B 



K*'mK* , 



rub 



{Kq |si(T^^6| B) = -'^{PtJiP'u-PuP^ 



Mm' 



ms + rriK* 



(8) 

(9) 

(10) 

(11) 
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The theoretical part of the correlator (5) can be calculated with the help of OPE at 
short distance. Up to operators having dimension 6, it is determined by the bare loop 
and power corrections from operators with d = 3 {qq), d = 4 (G^), d = 5 ml{qq) and 
d — 6 {qq)'^ ■ After calculating all power corrections, we obtain that (G^) and {qq)'^ give 
negligibly small correction and for this reason we do not present their explicit expressions. 

The spectral density for bare loop in Eq. (7) can be obtained by replacing the denom- 
inator l/ijP' — m^) of the quark propagators by — 27ri5(p^ — w?) in the initial Feynman 
integral. After standard calculation we get 



pt{p\p'\Q 



3 is the color factor, A 



{rrib (2s'A - A'u) - rus {2sM - A?i) } , (12) 



s — ml, A' 



s —rrig, u 



where ■ 
— Ass'. 

The integration region for the perturbative contribution in Eq. (7) is determined from 
the condition that the argument of all three 5 functions might vanish simultaneously, the 
bounds of s and s' are determined from the following inequality: 



-1 < 



liA - 2sA' 



< +1 



(13) 



The contributions of higher states is parametrized as the corresponding spectral density 
starting from s > sq and s' > s'q invoking quark hadron duality. 

We now present the results of calculations of dimension-3 and -5 operators: 



n 



(5) 



rr 

-ml {qq) 



1 



1 



mt 



ml 2ml + — 2mi,ms + 2Q^ 



(14) 
,(15) 



ml,r' ^ 



where r — 

In our calculations we do not calculate 0{as) corrections to the bare loop, and for 
consistency we also neglect 0{as) corrections in leptonic decay constants fx* and Jb- 

Substituting Eqs. (8), (14) and (15) into Eq. (7) and performing double Borel trans- 
formation in the variables and p'^, we get the following sum rule for the form factor 



rriB 



^^K* mi, {rnl/Ml+m^JMl) f 



2fK*mK* fsml 



I (27r) 



^ / dscis'pT(s,s',Q2)e-(^/^i'+^'/M|) 



+ e 



-(ml/M^+m^Mi) 



- {qq) + ml {qq) + 



Ml 3M| 



+ 



+ 



4M| 



+ 



2ml - 2mi,ms + 2ml + ^ 



12Mi2M| 

For completeness we also present sum rules for the form factors /+ and /_ 



(16) 



+ 



fBm,% fK*mK* 



e 



{2^y 
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{ml/M^+ml/M^) 



1 1 2 /- \ /3m^(m5 - 3m^(mf,-ms) 

- {qq) [mb - m,) - —m^ (qq) I - 



2Mf 



+ 



2M| 



2(mb — 2ms) 2(2m6 — m^) (mt — m5)(2m^ + mbrng + 2m^ + 2Q^ 



Ml 



+ 



1 ^m|/M2+m^,/M| 



+ 
+ 

where 



^ ^-{ml/M^+ml/M^) 

3m^(mb + ms) 2(mfe + Sm^) 2(3mb + m 



3ml{mb + rus) 



(17) 



(27r) 



^ (gg) (mfe + m,) + {qq) 



2Mt 



2M^ 



Ml 



'-2 ^"2 Mi 

{rUb + ms) {2ml + + '^'^l + 2Q^) ' 

m^m| 



P+ 



P- 



^{(A' + A)(l + A + 5) + [{mb + msf + g^](A + B)] 
^{[A' + A + K + m3)^ + g'](A-i?) + A'-A} , 



(18) 

(19) 
(20) 



and 



A = -[2s'A-mA'] , 
A 

B = i[2sA'-iiA] , 
A 

Using Eqs. (l)-(3), we get the following expression for the differential decay width: 



(21) 



dT{B K*l 
ds 



C^a^m^ 
30 727r5 



+ \CioU\ 



\VtbV;fv^X{l,ml.,s)< 



1 + rriKx 



(3 - ^;2)A(1, , §) + 12m\ (2 + 2m\. - s) |/+ 



+ 2(l-7n^.)Re[/+/:]+^|/_| 



C 
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(22) 



where 



4m| 



2 ' 

m% 



mb 



mb 
ruB 



me, 

me = , 

mB 



mK* = 



mx* 
mB 



and A(l, m|^., s) = 1 + m^-* + — 2s — 2m|-. (1 + s). 
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3 Numerical results of the sum rules 



In this section we present our numerical study for the form factors /+(?^), /-(?^), and fT{(f) 
for the B — > (1430)£'*"^~ decay. Prom the expressions of the sum rules for the form factors 
wc sec that the main input parameters of all form factors are quark condensates, leptonic 
decay constants fs and fx* of B and mesons, respectively, Borel parameters and 
M|, as well as continuum thresholds sq and Sq. 

The value of the quark condensate Bi ji—l GeV is {qq) = —(240 ± 10 MeV)^, mc{iji — 
rric) = (1.275 ± 0.015) GeV, = (4.7 ± 0.1) GeV [22]. The leptonic decay constant 
of B and the continuum thresholds Sq and Sq are determined from the analysis of two- 
point sum rules which have the values: ]b = 180 MeV [16], sq = (35 ± 2) GeV^ and 
Sq — (4.4 ± 0.4) GeV^ [23]. Few words about the value of the leptonic decay constant fx* 
are in order. This constant is calculated within the 2-point QCD sum rules method in [23] , 
including 0{as) corrections. As has already been noted, we neglect 0{as) corrections in 
the bare loop calculations, and hence, for consistency, we shall neglect them in calculation 
of the leptonic decay constant fx*- Neglecting 0{as) corrections from the result of [23], 
we obtain for the leptonic decay constant fx* = 340 MeV, which we shall use in further 
numerical calculations. 

The Borel parameters and M| are not physical quantities. The results for the phys- 
ically measurable quantities should be independent of them if the OPE can be performed 
up to infinite order. But, in sum rules, OPE is truncated to some finite order and for this 
reason the Borel parameters have to be chosen in such a working region that the physical 
results are practically independent on them. In choosing the working region of Mf and M| 
the conditions must be satisfied: 1) the contribution of the excited should be small, and, 
2) the power corrections should converge. 

Our numerical analysis shows that both conditions are satisfied and the best stabilities of 
all form factors are achieved when Mf and M| vary in the regions 8 GeV^ < Ml < 15 GeV^ 
and 2.5 GeV^ < M| < 4.5 GeV^, respectively 

The values of the form factors at = are 



where the errors are due to the variation of Borel parameters, the continuum thresholds 
So and s'q, the uncertainty in the condensate parameters, the variation of quark mass and 
meson decay constants. 

Here we would like to make the following cautionary note. Prom the sum rules expres- 
sions for the form factors it is seen that they are very sensitive to the magnitude of the 
leptonic decay constant fx* of the Kq{1A30) meson. The above-mentioned results for /i(0) 
are obtained at fx* = 340 MeV. But, if we use fx* = 300 MeV which is given in [14] 
(after normalizing it to our definition) we obtain 



/+(0) 
/-(O) 

MO) 



-0.31 ±0.07 



0.31 ±0.08 



0.26 ±0.07 



(23) 



/+(0) = 0.34 ±0.08, 
/_(0) = -0.34 ±0.07 , 
/y(0) = -0.29 ±0.07 . 



(24) 
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When compared with our results, /+(0) and /t(0) are in good agreement within the error 
limits, while the value of /-(O) is 50% larger compared to the prediction of [14]. 

In calculating the total width of the B ii'o (1430)£^£^ decay, wc need to know the 
dependence of the form factors f+{q'^), f-{(f) and friq'^) in the physical region 4m| < 
q^ < {niB — mK*Y- The q^ dependence of the form factors in the physical region can be 
calculated directly from sum rules, which has comprehensively been discussed in [19, 20] 
and [24]. We are restricted to consider a region of g^, calculated from the QCD side, where 
the correlator can be reliable. Our numerical analysis shows that this region is bounded 
as < < 8 GeV'^. In order to extend the results to full physical region, we look such a 
parametrization of the form factors that they coincide with the sum rule predictions in the 
above-mentioned region of q^. 

Our analysis shows that the best fit for the q^ dependence of the form factors can be 
written in the following form: 



1 — aiS + hiS^ 



(25) 



where i — -\-, — or T and s — q^ /m\. The values of the parameters /i(0), and hi are 
given in table 1. 





m 


ai 


h 


/+ 


0.31 ± 0.08 


0.81 


-0.21 


/- 


-0.31 ±0.07 


0.80 


-0.36 


h 


-0.26 ±0.07 


0.41 


-0.32 



Table 1: Form factors for B 
340 MeV. 



i^'Q(1430)£+^ decay in a three-parameter fit for fx* 



Using the equation of motion the form factor /t can be related to /_ as follows: 

(mb - ms) 



Mq' 



{ruB - rriK*) 



f- 



(26) 



We see that, within the errors Eq. (26) is in agreement with the computed form factor fx 
from QCD sum rules. 

The kinematical interval of the dilepton invariant mass q^ is 4m^ < Q"^ < (^s — ^k*)"^ 
in which the long distance effects (the charmonium resonances) can give substantial contri- 
bution. The dominant contribution to the B — > Kq{1A?,{))1'^1~ decays comes from the two 
low lying resonances J/V' and ■0', in the interval of 8 GeV"^ < < 14 GeV"^. In order to 
minimize the hadronic uncertainties we discard this subinterval by dividing the kinematical 
region of q^ as follows: 
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I 4m? < < (mj^ - 0.02 GeVf , 



II {mj^ + 0.02 GeVf < < {m^> - 0.02 GeVf , 

III (m^/ + 0.02 GeVf < < {rriB - mK*f . 

In Fig. (1) we present the dependence of the differential branching ratio for the B 
(1430)e"'"e~ decay on g^, as well as its dependence on q^ due solely to short distance 
effects (ae = case). 

Taking into account the q^ dependence of the form factors given in Eq. (25), performing 
integration over s, and using the total hfe time tb — 1.53 x 10~^^ s [25], we get the following 
for the branching ratios when only short distance contribution is taken into account: 

B{B X*(1430)e+e-) = 2.09 x 10"^ (2.68 x 10"^) , 
B{B X*(1430)/xV~) = 2.07 x 10"^ (2.66 x 10"^) , 
B{B K*{Um)T+T-) = 1.70 X lO-'^ (2.20 x lO"'^) , 

where the values in the parenthesis correspond to the choice of fx* = 300 MeV . 

Taking long distance effects into account in the above-mentioned kinematical regions, 
we get the following for the branching ratios: 



B{B ir*(1430)e+e-) = < 



1.97 X 10 region I , 

1.66 X 10~^ region II , at 33 = 1, 

3.90 X 10-1° region III , 



and 



' 2.13 X 10-'^ region I , 



B{B ir*(1430)e+e-) = < 



1.63 X 10-^ region II , at 33 = 2, 



3.38 X 10-1" region III , 



Xo*(1430)e+e- and 



at Jk* = 340 MeV . Note that the branching ratios of the B 
B ^'0(1430)//"'";^- decays are practically the same. 

It follows from these results that the dominant contribution comes from region I (low 
invariant mass region), and this can be attributed to the existence of the factor \/q^. Since 
at LHC-b lO^i-lO^^ pairs are hoped to be produced, the expected number of events for 
the B (1430)e"'"e- decay in the low invariant mass region is of the order of 10^-10^. 

Therefore, this mode sounds to be quite measurable at LHC-b. Comparing the value of 
the branching ratio for the B fCg (1430)e+e^ decay for the case when only short distance 
effects are taken into account, with the value calculated in region I, it can be said that 
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they practically coincide. Hence, in the light of above results, we can conclude that long 
distance effects do not give significant contribution in region I, and therefore, measurement 
of branching ratio in this region allows us to check the short distance structure of the 
effective Hamiltonian. 

The smallness of the value of B{B — > Kq{1A30)t'^t~) can be attributed to the small 
phase volume of this decay. Taking into account the efficiency for detecting r leptons, the 
measurement of the branching ratio for B A'q (1430)r+r~ is very difficult, even at LHC. 

As the concluding remark we can state that, we study the semileptonic rare B — > 
K^{U30)e+e- decay. The transition form factors of the B K^{U30)i+e- are calculated 
in the framework of three-point QCD sum rules. The branching ratios of the relevant 
decay for i = e, fi,T leptons, when short and long distance effects are taken into account, 
are estimated. From these results we conclude that the B Kq[14:30)£~^£~ {£ — e, //) decay 
can be measured in future planned experiments at LHC. 
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Figure captions 

Fig. (1) The dependence of the differential branching ratio for the B — > (1430)6+6" 
decay on g^, for the values of the fudge factor ae = (corresponding to short distance 
effects), ee — 1 and ee — 2. 
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Figure 1: 
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